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We present a theoretical analysis of magnetic toroidal moments in periodic systems, in the limit 
in which the toroidal moments are caused by a time and space reversal symmetry breaking arrange- 
ment of localized magnetic dipole moments. We summarize the basic definitions for finite systems 
and address the question of how to generalize these definitions to the bulk periodic case. We define 
the toroidization as the toroidal moment per unit cell volume, and we show that periodic boundary 
conditions lead to a multivaluedness of the toroidization, which suggests that only differences in 
toroidization are meaningful observable quantities. Our analysis bears strong analogy to the "mod- 
ern theory of electric polarization" in bulk periodic systems, but we also point out some important 
differences between the two cases. We then discuss the instructive example of a one-dimensional 
chain of magnetic moments, and we show how to properly calculate changes of the toroidization 
for this system. Finally, we evaluate and discuss the toroidization (in the local dipole limit) of four 
important example materials: BaNiF4, LiCoP04, GaFeOa, and BiFeOs. 



I. INTRODUCTION 

The recent resurgence of interest in magneto-electric 
multiferroics has prompted discussion of the relevance of 
the concept of magnetic toroidal moments in such sys- 
tems (see e.g. Refs. A magnetic toroidal 
moment is represented by a time-odd polar (or "axio- 
polar", see Ref. H) vector, which changes sign under 
both time inversion and space inversion, and is gener- 
ally associated with a "circular" or "ring-like" arrange- 
ment of spins (see Fig. [1] for some examples) Materials 
in which the toroidal moments are aligned cooperatively 
- so-called ferrotoroidics - have been proposed to com- 
plete the group of primary ferroics.-'-'- Further interest 
stems from the fact that the toroidal moment is related 
to the antisymmetric part of the linear magneto-electric 
tensor jiSiiii this points to an important role played by the 
magnetic toroidal moment for magneto-electric coupling 
phenomena. 

Magnetic toroidal moments in condensed matter sys- 
tems were first studied in the former Soviet Union during 
the 1980s, mostly in the context of the so-called "exci- 
tonic insulator" model (see Refs. lolfl^ for reviews of this 
work). At about the same time, Sannikov and Zhelu- 
dov proposed the toroidal moment as the primary or- 
der parameter for the low temperature phase transition 
in multiferroic nickel iodine boracite.^'^ Since in these 
early studies the toroidal moment was mostly treated as 
a macroscopic order parameter, no particular attention 




FIG. 1: Simple arrangements of magnetic moments which can 
lead to toroidal moments, (a) and (b) have equal and opposite 
toroidal moments. The antiferromagnetic arrangement in (c) 
has a toroidal moment, whereas that in (d) does not. 



was paid to the peculiarities arising from the microscopic 
definition of the toroidal moment. It is therefore the pur- 
pose of the present paper to give a detailed analysis of the 
properties of toroidal moments starting from the micro- 
scopic definition and focusing especially on effects result- 
ing from the periodic boundary condition in crystalline 
solids. In particular, we address the following questions: 

1. How should the toroidal moment density, or 
toroidization, of a bulk periodic solid be formally 
defined? 

2. Is there a consistent way to treat the origin depen- 
dence of the toroidal moment? 

3. What are the consequences of the periodic bound- 
ary conditions within a bulk crystalline solid? 

In addition, we apply our newly developed concepts to 
evaluate and analyze the toroidization of four example 
materials: the antiferromagnetic ferroelectrics BaNiF4 
and BiFeOa, the polar ferrimagnet GaFeOs, and the 
strongly magneto-electric material LiCoP04. 

To avoid confusion, we point out that there also has 
been some recent discussion about electric toroidal mo- 
ments g, defined as g = ^ J^i '"'i ^ Pi where pi is the local 
electric dipole moment at position and the summa- 
tion extends over all dipole moments in the system.— d^d^ 
The vector g is fundamentally different from the mag- 
netic toroidal moment, since it is both time- and space- 
inversion symmetric. It has been used to characterize cir- 
cular domains in nano-scale ferroelectric materials f^^i^^ 
In the following we exclusively discuss the case of mag- 
netic toroidal moments, i.e. the term "toroidal moment" 
is always used in the sense of "magnetic toroidal mo- 
ment". We point out, however, that some of our gen- 
eral considerations regarding origin-dependence and the 
effect of periodic boundary conditions are applicable to 
the case of electric toroidal moments as well. 

This paper is organized as follows. We begin by sum- 
marizing some of the basic definitions and then discuss 
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the hmit where the toroidal moment is caused by a time 
and space reversal symmetry breaking arrangement of lo- 
calized magnetic moments (Sec. HIl. In Sec. IIIII we then 
analyze the origin dependence of the toroidal moment 
by decomposing the magnetic moment distribution into 
a fully compensated, generally non-coUincar, antiferro- 
magnetic part and a non-compensated, coUinear, ferro- 
magnetic part. In Sec. IIVI we define the toroidization 
as the toroidal moment per unit cell volume, and we 
show that the periodic boundary conditions lead to a 
multivaluedness of the toroidization, similar to the case 
of the electric polarization in bulk periodic solids (see 
Refs. This multivaluedness suggests that, as 

in the case of the electric polarization, only differences in 
toroidization can be physically observable quantities. In 
SeclVlwe illustrate some consequences of the periodicity 
by using the example of a one-dimensional antiferromag- 
netic chain of magnetic moments. In Sec. IVII we evaluate 
the toroidizations of four example materials: BaNiF4, 
LiCoP04, GaFeOa, and BiFeOg. Finally, in Scc.|VlIlwe 
summarize our main conclusions and discuss the corre- 
spondence between the microscopic toroidal moment de- 
scribed in this paper and some phenomenological quanti- 
ties with the same symmetry, that have recently appeared 
in the multiferroics literature. 



II. DEFINITIONS 

The toroidal moment t corresponding to a current den- 
sity distribution j{r) is defined as (see Ref. 0): 

(^(^•j)-2r2j)dV , (1) 

where c indicates the speed of light in vacuum. The 
toroidal moment in the form of Eq. ([T]) emerges from 
the multipole expansion of an arbitrary localized current 
distribution.^ Its physical significance can be seen by not- 
ing that Eq. ^ is identically satisfied by the current dis- 
tribution 

j{r) = cV X V X 5{r)t , (2) 

which represents an elemental toroidal moment centered 
at the origin. Evaluating the interaction energy E of this 
current density with the electromagnetic field A(r) yields 
(after partial integration): 

E = -^J j{r) ■ A{r) d\ = -t ■ V x B{0) , (3) 

where B{0) — V x A{r)\r=Q is the magnetic field at the 
site of the toroidal moment. From Eq. ([3]) it can be seen 
that the toroidal moment couples to the curl of the mag- 
netic field such that a toroidal system in a magnetic field 
has lowest energy when its toroidal moment is aligned 
parallel to the curl of the magnetic field. 

The definition of the toroidal moment can be re-cast 
into a more convenient form, by noting that the current 



vector can be decomposed into longitudinal (V x j|| = 0) 
and transversal {V ■ j± — 0) parts. The longitudi- 
nal part of j{r) is related to time derivatives of the 
charge multipole moments through the continuity equa- 
tion p+V -j = 0, and does not contribute to the toroidal 
moment)^ The transverse part of the current density 
j±{r) can be written as the curl of the magnetization 
density /.t(r)ii^ 

j±{r) = cV X fi{r) . (4) 

Inserting Eq. ^ in Eq. ^ gives the toroidal moment in 
terms of /i,(r):^" 

* = i y r X fi{r) d^r . (5) 

While in principle, for a finite system with known magne- 
tization density, this expression can be used to evaluate 
the toroidal moment, it is not directly applicable to ex- 
tended systems where periodic boundary conditions are 
employed. The difficulties resemble those encountered in 
early attempts to calculate the electric polarization (see 
Ref . fisl) : for a general continuous magnetization density 
/x(r), Eq. (O evaluated over one unit cell will lead to 
arbitrary values, depending on the special choice of unit 
cell used in the calculation. 

In the case of the electric polarization, a general so- 
lution to this problem is achieved by evaluating the 
electric polarization directly from the electronic wave- 
function using Wannier representations] ""^^i^^ i^^ In prin- 
ciple, a similar approach seems appropriate for the case 
of the toroidal moment. In the present paper, we pur- 
sue a somewhat simpler yet rather instructive approach 
by assuming that the magnetization density can be well 
represented by a distribution of localized magnetic mo- 
ments {nia} at sites r^: 

Hiocir) ^^maS{r ~ ra) ■ (6) 

a 

This results in the toroidal moment: 

t = i ^ X ma . (7) 

a 

The simplification to localized magnetic moments avoids 
the technical difficulties of a full gauge invariant wave- 
function formulation but retains all the peculiarities re- 
sulting from the periodic boundary conditions within 
bulk systems. It also allows a more intuitive analysis 
of several important prototype systems. Our study thus 
represents a first step towards a full microscopic theory 
of toroidal moments in crystalline solids and can be used 
as the basis for future developments. 

We note however, that in some cases the restriction 
to localized magnetic moments can represent a severe 
simplification. Obviously, the local moment picture ne- 
glects the possibility of toroidal moments arising from 
non-localized magnetization densities, but in addition 
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FIG. 2: Decomposition of a ferrimagnetic arrangement of two 
localized moments (left) into its fully compensated component 
(middle) and its uncompensated "ferromagnetic" component 
(right). The ferromagnetic component at each site is the total 
magnetic moment divided by the total number of moments, 
rha = \{m\ + m2), and the compensated part is the differ- 
ence between the magnitude of the original local moment and 
the uncompensated contribution. 



the symmetry of the magnetic moment distribution in 
Eq. ([6|) can eventually be higher than the full magnetic 
space group symmetry represented by the original mag- 
netization density fJ-^r). This can occur even in systems 
that are usually well described in terms of localized mag- 
netic moments (see for example our discussion of BiFeOa 
in Sec. IVI D|) . The localized moment approach also ne- 
glects the possibility that the localized current distribu- 
tion around the atomic site, which gives rise to the lo- 
cal magnetic dipole moment, simultaneously gives rise 
to a localized toroidal dipole moment. Such "atomic" 
toroidal moments have been discussed in the context of 
atomic multipole moments and can in principle be mea- 
sured by resonant x-ray spectroscopy.— Here, we restrict 
our discussion to the case of toroidal moments "on the 
unit cell scale" and disregard the possibility of toroidal 
contributions "on the atomic scale" . 

Using Eq. ([7]) we can straightforwardly evaluate the 
toroidal moments of the arrangements shown in Fig. [1] 
Taking the horizontal magnetic moments to be spaced a 
distance a apart along the y direction, and the vertical 
moments a distance a apart along x, the toroidal mo- 
ments of arrangements a) and b) in Fig. [1] are t = ±asz, 
where s is the magnitude of the individual magnetic mo- 
ments. The toroidal moment of Fig. [it is t = —^asz, 
whereas it is zero for the arrangement shown in Fig. [TJl, 
since in this case the moment vectors are aligned parallel 
to the vector connecting the two sites. 



III. ORIGIN DEPENDENCE 

It can easily be seen that for systems with non- 
vanishing magnetic dipole moment, 

m= ^ J r X j{r)d^r = J fj,{r)(fr =::^ma , (8) 

a 

the toroidal moment in Eqs. ((T|), (O, and ^ depends on 
the choice of origin. In particular, for a change of origin 
defined by r ^ r' = r + i?o the toroidal moment changes 
a.st^t' = t + l/2Ro X m. 



To further analyze this origin dependence, we decom- 
pose the original magnetic moment distribution {m^} 

into two parts: a totally compensated part {mi"'' = 

nia — m/N}, with no net magnetization, J2a ''^'a^ = 0, 
and an uncompensated "ferromagnetic" part {rha — 
m/N}, where m — J2a''^a is the total magnetic mo- 
ment and TV is the total number of localized moments 
(see Fig. [5] for a simple example). This results in a cor- 
responding decomposition of the toroidal moment into: 

a 

and 

i=^Rxm , (10) 

where R = 1/-/V^^ is the average magnetic moment 
position. By construction, only the toroidal moment t 
depends on the choice of origin, whereas the part t^"-* is 
origin independent. 

For a macroscopic toroidal moment to occur, the mag- 
netic moment distribution has to break both time and 
space reversal symmetries. In the compensated moment 
distribution {rria^} this can happen in several ways, de- 
pending both on how the magnetic moments are ori- 
ented and on how they are positioned. On the other 
hand the uncompensated distribution {rha} provides less 
freedom. Due to the non-vanishing magnetic dipole mo- 
ment, the configuration {rha} always breaks time rever- 
sal symmetry. However, the only possibility for such a 
"ferromagnetic" moment distribution to simultaneously 
break space inversion symmetry is that the magnetic mo- 
ments are positioned in a non-centrosymmetric way. A 
nonzero toroidal moment of the uncompensated part of 
any moment distribution is therefore always related to an 
inversion symmetry-breaking arrangement of the under- 
lying ionic lattice, whereas this does not necessarily have 
to be the case for the compensated moment distribution 
{rria^}, where the inversion symmetry can also be lifted 
by the orientation of the various magnetic moments. -^^ 

It will become clear in the following section, that only 
differences in toroidal moment should have any physi- 
cal significance. In the case of the origin-dependent part 
t, such differences in toroidal moment must be related 
to corresponding displacements of the moment positions, 
and the change in the toroidal moment is then given by 
A* = 1/2 AR X m, with = l/NJ2a and Ar^ 
being the displacements of the individual magnetic mo- 
ments. Thus, if a consistent choice of origin is used for the 
initial and final configuration, the corresponding change 
in the toroidal moment is a well-defined physical quan- 
tity. For example, if the initial reference configuration 
is centrosymmetric, then the change in toroidal moment 
resulting from a symmetry-breaking structural distortion 
can be interpreted as the spontaneous toroidal moment 
of the system. 
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Earlier applications of Eq. ([7]) did not perform the ex- 
plicit decomposition described above, but instead eval- 
uated the toroidal moment with respect to the "center 
of the unit cell",^"^ without specifying exactly how this 
center of the unit cell is defined. We point out that the 
origin dependent contribution to the toroidal moment t 
vanishes, if the "center of the unit cell" defined by R is 
taken as the origin,^'* However, it is important to realize 
that for cases where the toroidal moment changes as a re- 
sult of a structural distortion, R in general also changes 
(see the discussion in the previous paragraph). In such 
cases the origin should be taken to be the same for both 
structural modifications. 

Finally, we note that in the case of a non- localized mag- 
netization density, the "uncompensated" part of /i(r) 
would correspond to a uniform, i.e. perfectly homoge- 
neous (r-independent), magnetization density p, — m/V, 
where V is the total volume of the system. Since such 
a perfectly homogeneous magnetization density fi can 
never break space-reversal symmetry, it does not con- 
tribute at all to the toroidal moment of the system. 
Thus, in a non-localized description, the decomposition 
into compensated and uncompensated parts results in a 
perfect separation between dipolar and toroidal contri- 
butions to the magnetization density. Due to the "in- 
homogeneity" that is intrinsic to the localized moment 
description, the decomposition is not fully complete in 
this case, and the toroidal contribution of the "ferromag- 
netic" part has to be analyzed separately. 

In our explicit examples in Sees . IVl and I VII we will only 
consider systems with fully compensated moment con- 
figurations and postpone the further analysis of toroidal 
moments resulting from "ferromagnetic" configurations 
to future work. 



IV. THE CASE OF A PERIODIC BULK 
SYSTEM 

We now proceed to the case of bulk periodic crystals 
with an infinite periodic arrangement of magnetic mo- 
ments. We begin by defining the toroidal moment per 
unit volume or toroidization T = t/V, where V is the 
volume of the system with toroidal moment t. Then, 
for a large finite system containing N identical unit cells 
each of volume f2: 



2NQ. 



Eiri + Rn) X mi 



2Nn 



(11) 

(12) 



Here, are the positions of the magnetic moments rrii 
relative to the same (arbitrary) point within each unit 
cell, Rn is a "lattice vector" with index n, and we have 
used the fact that the orientation of the magnetic mo- 
ments is the same in each unit cell. The summation over 
i indicates the summation over all moments within a unit 



cell, and that over n indicates the summation over all unit 
cells. Expanding the cross product, we obtain: 



, . , X rrii -] 



^^i?„xm (14) 



^^r, xm, , (15) 



2n 



where in the last step we have assumed that the sum 
over all lattice vectors contains both i?„ and —Rm so 
that Rn — 0. This is true for any infinite Bravais 
lattice. Thus, the toroidal moment of a system of N unit 
cells is just N times the toroidal moment evaluated for 
one unit cell, and the corresponding toroidizations are 
identical. 

In an infinite periodic solid, we have a freedom in 
choosing the basis corresponding to the primitive unit 
cell of the crystal. In particular, we can translate any 
spin of the basis by a lattice vector Rn without changing 
the overall periodic arrangement. However, such a trans- 
lation of magnetic moment rrii by i?„ leads to a change 
in the toroidization as follows: 



(16) 



The freedom in choosing the basis corresponding to the 
primitive unit cell thus leads to a multivaluedness of the 
toroidization with respect to certain "increments" (de- 
fined by Eq. ^TE\i ) for each magnetic sub-lattice i and 
lattice vector i?„. 

This multivaluedness of the toroidization is in 
strong analogy to the modern theory of electric 
polarization,iS]i2ii^ where the polarization changes by 
ei?„/r2 if one translates an elementary charge e by a 
multiple of a lattice vector Rn- The resulting multi- 
valuedness has led to the concept of the "polarization 
lattice" corresponding to a bulk periodic solid,— with 
ei?„/ri called the "polarization quantum" if Rn is one of 
the three primitive lattice vectors. Eq. suggests the 
existence of an analogous "toroidization lattice". How- 
ever, the vector product in Eq. and our classical 
treatment of magnetic moments can lead to arbitrary 
projections of rrii on a certain axis, and therefore the 
structure of the "toroidization lattice" is quite different 
from that of the polarization lattice. In the most general 
case, if there are r magnetic basis atoms in the primitive 
unit cell, there can be 3r linearly independent toroidiza- 
tion increments. This can lead to cases, where multiple 
incommensurate increments exist along certain crystal- 
lographic directions. In addition, for a coUinear moment 
configuration no toroidization component is allowed par- 
allel to the global magnetic axis. Thus, the set of al- 
lowed toroidization values does not necessarily have the 
same translation symmetry as the corresponding crystal 
structure and does not necessarily form a Bravais lattice. 
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whereas this is always true in the case of the electric 
polarization. In practice, the magnetic symmetry of the 
system can significantly reduce the number of linearly in- 
dependent toroidization increments (see Sec. IVII for some 
realistic examples. 

In spite of the difficulties associated with the multi- 
valuedness of the polarization, it is now recognized that 
only differences in the polarization lattices between dif- 
ferent ionic configurations are in fact measurable quanti- 
ties, such as for example the difference between two oppo- 
sitely polarized states of a ferroelectric crystal or between 
a centrosymmetric non-polar reference structure and the 
actual polar crystal. These differences are the same for 
each point of the polarization lattice and are thus well- 
defined quantities. In analogy with the case of the electric 
polarization we suggest that only differences in the set of 
allowed toroidization values, corresponding to two differ- 
ent bulk configurations, are physically observable quanti- 
ties, such as the difference in toroidizations between two 
domain states of a ferrotoroidic, or the difference between 
a ferrotoroidic state and its non-toroidic paraphase. Such 
quantities can be obtained by monitoring the change in 
toroidal moment on one arbitrarily chosen branch within 
the allowed set of values, when transforming the system 
from the initial to the final state along a well-defined 
path. 

Finally, we emphasize that the multivaluedness of the 
toroidization and the possible origin dependence of the 
toroidization are two independent features with different 
origins. Both features are rooted in the fundamental def- 
inition of the toroidal moment in terms of the position 
operator r, but whereas the origin dependence appears in 
both finite and infinite systems if there is a non- vanishing 
magnetic dipole moment, the multivaluedness is caused 
by the periodic boundary conditions in a bulk solid, and 
is independent of an eventually non-vanishing magnetic 
dipole moment. 



V. A ONE-DIMENSIONAL EXAMPLE 

A. The periodic non-toroidal state 

To illustrate some consequences of the multivalued- 
ness of the toroidization in periodic systems described 
in the previous section, we now consider the example 
of a one-dimensional antiferromagnetic chain of equally 
spaced magnetic moments as shown in Fig. [3^. The mo- 
ments, with magnitude s, are spaced a distance a apart 
from each other along the y axis, and are alternating in 
orientation along ix. Thus, the unit cell length is 2a and 
there are two oppositely oriented magnetic moments in 
each unit cell. Since this configuration does not possess a 
macroscopic magnetic dipole moment, the corresponding 
toroidal moment is origin-independent, and a decompo- 
sition into compensated and uncompensated parts is not 
required. 

The arrangement of magnetic moments in Fig. [3^ is 



a) *i 



b) 



FIG. 3: Calculation of the toroidization for two different one- 
dimensional antiferromagnetic periodic arrangements of mag- 
netic moments. Our choice of unit cell is indicated by the 
shaded area in each case, a) shows a non-toroidal state, which 
is space-inversion symmetric with respect to each moment 
site, b) is a toroidal state. 



space-inversion symmetric with respect to each moment 
site and thus cannot exhibit a macroscopic toroidal mo- 
ment. Furthermore, even though the arrangement in 
Fig. [3^ breaks time reversal symmetry microscopically^ 
there exists a symmetry transformation which combines 
time inversion with a translation of all moments by the 
distance a along the y direction. According to Neumann's 
principle (see e.g. Ref. [25h . the macroscopic properties 
of a system cannot depend on such microscopic transla- 
tions, i.e the macroscopic properties are determined by 
the point group of the system and not by its space group. 
Therefore, time reversal symmetry is not broken macro- 
scopically for the moment configuration in Fig. [3^ and 
its point group contains time inversion as a symmetry 
element. No macroscopic toroidal moment can thus re- 
sult from this configuration, in spite of the fact that an 
isolated unit cell would exhibit a toroidal moment. 

The toroidal moment of the single unit cell highlighted 
in Fig. [3^, calculated using Eq. ([7]) , is identical to that 
calculated for the finite spin configuration in Fig. [!};, 
i.e. t = —^saz, and the corresponding toroidization, 
T = t/Vt — —jZ (since the "volume" ft of the one- 
dimensional unit cell is just its length, 2a). The elemen- 
tary toroidization increment in this case is AT = ±|i, 
which means that the toroidization of the unit cell is ex- 
actly equal to one half of the toroidization increment, and 
the allowed toroidization values for the periodic arrange- 
ment are T„ = + n)^z, where n can be any integer 
number. 

We see that in our example the allowed toroidiza- 
tion values form a one-dimensional lattice of values, cen- 
trosymmetric around the origin. This is analogous to the 
case of the electric polarization, where the polarization 
lattice is invariant under all symmetry transformations of 
the underlying crystal structure. In particular, the polar- 
ization lattice corresponding to a centrosymmetric crys- 
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tal structure has to be inversion symmetric. This can be 
achieved either by a lattice that includes the point P — 0, 
or by the same lattice but shifted from the origin by ex- 
actly half a polarization quantum. We see that the same 
holds true for the toroidization of our one-dimensional 
example, so that a centrosymmetric set of toroidization 
values can be understood as representing a non-toroidal 
state of the corresponding system. 

We emphasize again, however, that the toroidization 
case is slightly more involved than the case of the electric 
polarization. In particular, due to the fact that several 
incommensurate toroidization increments can exist along 
the same crystallographic direction, there are many more 
allowed values of the toroidization in the non-toroidal 
state than for the polarization in the non-polar state, as 
can be seen for the system BaNiF4 discussed in Sec. I VI Al 

B. Toroidal state and changes in toroidization 

In order to obtain a nontrivial "macroscopic" 
toroidization the system has to break both space and 
time inversion symmetry. In the case of the one- 
dimensional antiferromagnetic chain this can be achieved 
by "spin pairing" , i.e. if the distances between neigh- 
boring magnetic moments alternate as shown in Fig. ^jp. 
Here the magnetic moments of magnitude s are spaced 
alternately a distance of (1 — d)a and (1 + d)a apart from 
each other along the y axis (— 1 < d < 1), and again are 
alternating in orientation along ix. The non-toroidal 
example above corresponds to d = 0. Since the unit 
cell size is the same as in the non-toroidal case, the el- 
ementary toroidization increment is again AT = ±|£. 
The toroidization of the unit cell indicated in Fig. [Sb is 
T = — (1 — d)|z, so that the allowed values of T for the 
full periodic arrangement are: 

T,M)={n-^^y^z . (17) 

Fig.[3]shows the allowed toroidization values as a function 
of the displacement d of the spins from their positions in 
the centrosymmetric, non-toroidal state. 

The change in toroidization between two configura- 
tions with d — di and d = d2 for a certain branch n 
is given by: 

T„(d2)-r„(di)-|(d2-di)i , (18) 

i.e. it is independent of the branch index n. In partic- 
ular, if the non-centrosymmetric distortion is inverted 
(di = do, d2 = —do), the change in toroidization is 
2Ts = ^ so that Ts = ^ can be interpreted as the spon- 
taneous toroidization^ again in analogy to the case of the 
electric polarization, where the spontaneous polarization 
is given by the branch-independent change in polariza- 
tion compared to a centrosymmetric reference structure. 

Another possible way to alter the toroidization is by 
changing the orientation of the magnetic moments in- 
stead of changing their positions. In particular, we expect 



T(sz) 




FIG. 4: Allowed values of the toroidization for the antiferro- 
magnetic chain of Fig.|3]as a function of displacement d from 
the non-toroidal case (d = 0). The cartoons at the bottom in- 
dicate the corresponding positions of the magnetic moments 
within the unit cell. 

that a full 180° rotation of all magnetic moments, which 
is equivalent to the operation of time reversal, should 
invert the macroscopic "spontaneous toroidization" , and 
should therefore lead to the same change 2Ts as discussed 
above. If we allow the magnetic moments to rotate out 
of the X direction, while preserving the antiparallel align- 
ment of the two basis moments, the toroidization along 
the z direction is given by 

r„^(d,a)= (^n-i^)|cosa , (19) 

where a is the angle between the magnetic moments and 
the X direction. The change in toroidization for a full 
180° rotation of the moments is thus: 

T„^(do,180°)-r,:(do,0°) = -^ + s(2n+l) , (20) 

and apparently depends on the branch index n. However, 
if one calculates the same change in toroidization for the 
non-toroidal state with d = 0, one obtains: 

y„-(0,180°)-r,:(0,0°) = s(2n+l) . (21) 

Obviously, in this case the corresponding change in 
macroscopic toroidization should be zero, since both the 
initial and final states (and all intermediate states) cor- 
respond to a non-toroidal configuration and thus Tg = 0. 
If one subtracts the improper change in T^, Eq. (j2ip . 
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FIG. 5: Effect on the magnetic moment configuration of 
Fig. (center) of a reversal of all magnetic moments (right) 
and of a reversal of the non-centrosymmetric distortion d 
(left). Note that the right and left final states are identi- 
cal, with the moments on the left translated by half a unit 
cell compared to those on the right. 

from the change in toroidization calculated in Eq. (PO)) . 
one obtains the proper change in toroidization 2Ts = ^ , 
which is identical to one obtained by inverting the non- 
centrosymmetric distortion d. Here, we use the terminol- 
ogy "proper" and "improper" in analogy to the case of 
the proper and improper piezoelectric response,"^ where 
a similar branch dependence is caused by volume changes 
of the unit cell, and the improper piezoelectric response 
has to be subtracted appropriately. 

Fig.[5]shows the initial and final states for the two cases 
where either the moment displacements or the magnetic 
moment directions are inverted. The two final states are 
equivalent except for a translation of all moments by half 
a unit cell along y, which, again due to Neumann's prin- 
ciple, is irrelevant for the macroscopic properties. The 
spontaneous toroidization of the state on the left in Fig. [5] 
is therefore the same as for the state on the right side of 
Fig. [51 

VI. TOROIDIZATIONS FOR SOME EXAMPLE 
MATERIALS 

To further illustrate the concept of toroidal moments 
in crystalline solids and to investigate the consequences 
of the definitions and simplifications outlined in the 
preceding sections for real systems, we now evaluate 
the toroidizations of four example materials: BaNiF4, 
LiCoP04, GaFeOa, and BiFeOs. All these materials have 
been discussed recently in the context of multiferroics, 
magneto-electric coupling, or ferrotoroidics. BaNiF4 and 
BiFeOa are both antiferromagnetic ferroelectrics, with 
additional weak ferromagnetism in the case of BiFeOs, 
and the possible coupling between the various order pa- 
rameters in these systems has recently been studied us- 
ing first principles techniques i^s^S LiCoP04, which is 



TABLE L Positions rt and moment directions rrii of the mag- 
netic Ni cations in BaNiF4. a, b and c are the orthorhombic 
lattice constants, S represents an internal structural parame- 
ter, 9 is the angle between the magnetic moments and the or- 
thorhombic 6 axis, and S is the magnitude of the Ni magnetic 
moment. The corresponding lattice vectors are: ai = (a, 6,0), 



a2 = 


= (a/2. 


-i)/2,0), and 
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(0,0, c). 
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not ferroelectric, exhibits a rather large linear magneto- 
electric effect)^ and the observation of ferrotoroidic do- 
mains in this material using nonlinear optical techniques 
has been reported.^ GaFeOs is a magnetic piezoelec- 
tric which exhibits a strong asymmetry in the magneto- 
electric tensor. This asymmetry has been interpreted 
to result from a non-vanishing toroidal moment in this 
material.— 



A. BaNiF4 

BaNiF4 belongs to an isostructural family of antiferro- 
magnetic ferroelectrics with composition BaMF4, where 
M can be Mn, Fe, Co, or Nii^i^ ft was recently shown 
that BaNiF4 exhibits two distinct antiferromagnetic or- 
der parameters and that the order parameter correspond- 
ing to the "weak" (secondary) antiferromagnetic order 
can be reversed by using an electric field. The crystallo- 
graphic structure of BaNiF4 is orthorhombic correspond- 
ing to space group Cmc2i The magnetic space group 
is Pa2i,— which contains a non-primitive translation 
along the a direction combined with time inversioni^ 
This reflects the fact that the antiferromagnetic ordering 
in BaNiF4 leads to a unit cell doubling compared to the 
paramagnetic phase. The corresponding (macroscopic) 
magnetic point group is thus 21', which does not break 
time inversion, similar to the case of the one-dimensional 
antiferromagnetic chain discussed in Sec. lV Al Therefore, 
neither a macroscopic magnetization nor a macroscopic 
toroidal moment is allowed in this symmetry. We note 
that, in general, when the magnetic ordering leads to a 
unit cell doubling compared to the paramagnetic phase, 
then the system is always macroscopically time reversal 
symmetric and therefore non-toroidal. 

Nevertheless it is instructive to examine the effect of 
the periodic boundary conditions and the resulting mul- 
tivaluedness of the toroidization for this trivial case. The 
magnetic unit cell of BaNiF4 contains 4 magnetic Ni 
ions, whose positions and spin directions are listed in 
Table HI A decomposition into fully compensated and 
uncompensated components is not necessary, since there 
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is no macroscopic magnetization in this system. Applica- 
tion of Eqs. ([7]) and (|16p leads to a "toroidization lattice" 
of the form: 



COS0 



sm ( 



/ 




|(/ + m-2n) 
{l + ni + 2n- 




4) 



(22) 



where 9 is the canting angle of the magnetic moments (of 
magnitude S) relative to the orthorhombic b direction, 
and a, 6, c are the orthorhombic lattice parameters, k, 
I, m, and n are arbitrary integer numbers corresponding 
to four different independent toroidization increments in 
this system. Due to the pairwise coUinear spin structure 
in BaNiF4 the original 4 x 3 = 12 increments according 
to Eq. p6| are reduced to 2 x 3 = 6. Additional sym- 
metries reduce the number of independent increments to 
4. Due to the base-centered orthorhombic Bravais lattice 
of BaNiF4, these toroidization increments are in general 
neither parallel to the cartesian coordinate directions nor 
perpendicular to each other. It can be seen from Eq. (|^^ 
that for a general value of the angle 9, the allowed values 
of T do not form a Bravais lattice, and the multivalued- 
ness of T is more complex than in the case of the electric 
polarization. Nevertheless, the set of allowed toroidiza- 
tion values in BaNiF4 is inversion symmetric as required 
by the time-symmetric point group. In the absence of 
canting (9 — 0°), the toroidization and the toroidization 
increment along the b direction would be zero; however 
we see that for small 9 there is a small non-zero incre- 
ment of the toroidization along this direction, reflecting 
the small component of the magnetic moments perpen- 
dicular to b. Thus, the example of BaNiF4 shows that the 
structure of the set of allowed toroidization values is in 
general more complex than the "polarization lattice" in 
the modern theory of electric polarization in crystalline 
solids. 



B. LiCoP04 

The observation of ferrotoroidic domains in LiCoP04 
usingnonlinear optical techniques has been reported in 
Ref. 0. LiCoP04 crystallizes in the olivine structure 
with the orthorhombic space group Pnmaf^ and orig- 
inally it was believed that the magnetic moments of the 
four Co ions in the unit cell are antiferromagnetically 
aligned along the orthorhombic b direction. '^^ This mo- 
ment configuration corresponds to the magnetic space 
group Pnma' , which contains the operation of simulta- 
neous time and space inversion, but not the pure time 
and space inversion separately, and thus allows the ex- 
istence of a macroscopic toroidal moment. Furthermore, 
the toroidal moment is required to be aligned parallel 



TABLE II: Positions and magnetic moment directions rtii 
of the Co cations in LiCoP04, according to Ref. [S^ and as- 
suming a moment rotation towards the c direction, e and 5 
are internal structural parameters corresponding to Wyckoff 
positions 4c of the Pnma space group. 6 is the angle between 
the magnetic moments and the b axis and S is the magni- 
tude of the Co magnetic moment. Experimental values are 
e = 0.0286, S = 0.0207, and 6 ^ 4.6° M The orthorhombic lat- 
tice parameters are a = 10.20 A, 6 = 5.92 A, and c = 4.70 AM 
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to the c axis. Recently, it was found that the magnetic 
moments in LiCoP04 are rotated slightly away from the 
b direction by an angle 9 k, 4.6°, while preserving the 
overall coUinear magnetic structure.^- For a moment ro- 
tation within the b-c plane, the rotated moment config- 
uration corresponds to a lower symmetry with the mag- 
netic point group 2'/to, which allows a toroidal moment 
also along the orthorhombic b direction. At the moment, 
it is not clear what the primary order parameter for this 
additional symmetry lowering is, although since it does 
not change the relative antiferromagnetic arrangement of 
the spins a toroidal origin has been proposed.'' Further- 
more, a weak magnetization has been measured along the 
b direction^i^ which indicates an even lower point group 
symmetry of 2'. Since the small magnetization occurs 
along the h direction, i.e. parallel to the direction of 
the antiferromagnetic alignment, it has been described 
as weak ferrimagnetism. 

Here we calculate the toroidization corresponding to 
the fully compensated antiferromagnetic configuration 
with magnetic point group 2'/r7i, where all spins are ro- 
tated from the b direction towards the c direction by 
an angle 9. The corresponding positions and spin di- 
rections of the four magnetic Co cations within the unit 
cell are listed in Table [TTl The toroidization resulting 
from Eqs. ([7]) and (fT6|) is: 



cos ( 



(4eT0a ] 



sm ( 



bn 
(4e±0a 

V 







(23) 

Here, a, b and c are the orthorhombic lattice constants 
and n, to, and I are arbitrary integer numbers. It can 
be seen that there is only a trivial component of T along 
the a direction and that there is a nontrivial part of the 
toroidal moment proportional to e, which rotates from 
the c direction for 9 = Q° towards the b direction for 
9 ^ 90°. For e = the system is centrosymmetric 
and thus non-toroidal. The spontaneous toroidization 
is given by = ^fr- experimentally reported 

parameters listed in Table [TT] and using the formal mag- 
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netic moment S = 3/xb of Co^+ , the corresponding value 
is Ts = 6.17- 10~'^/iB/A (icos^ + ysin6'), correspond- 
ing to a spontaneous toroidal moment per unit cell of 
tg — 1.75/xbA (i cos + y sin6'). Thus, the magnetic mo- 
ment rotation described by 9 results in a corresponding 
rigid rotation of the toroidization. 



C. GaFeOg 

Ga2_2:Fea;03 was the first material that was found 
to exhibit both piezoelectricity and a macroscopic 
magnetization^ It was also the first known material 
with a spontaneous magnetization that simultaneously 
exhibits a linear magneto-electric effect.— The system 
Ga2_a;Fe^03 has been studied recently because of its in- 
teresting optical properties, which result from the simul- 
taneous breaking of both space and time inversion sym- 
metries in this material42ii^ In addition, an asymmetry 
of the magneto-electric tensor has been measured and 
was attributed to the existence of a toroidal moment in 
this systemj22, 

The crystal structure (space group Pc2in) of 
Ga2-xFex03 contains four inequivalent cation sites, one 
with tetrahedral oxygen coordination and three different 
octahedrally coordinated sites4^ The Fe cations occupy 
predominantly two of the octahedral sites (called Fel and 
Fe2), but there is also a sizable Fe occupation on the third 
octahedral site (Ga2), whereas the tetrahedral site (Gal) 
is occupied mainly by Ga. The exact occupation of the 
various cation sites depends on the composition x as well 
as on the preparation technique and sample history4^ 

Two different magnetic structures have been proposed 
for this system, both of which are consistent with the 
magnetic space group which allows for a macro- 

scopic magnetization along the crystallographic c direc- 
tion. Abrahams and Reddy originally proposed a canted 
antiferromagnetic structure with compensating magnetic 
moments within the a-b plane and a net magnetization 
along c4i On the other hand Arima et al. (Ref. |43| ) 
recently interpreted their data in terms of a coUinear fer- 
rimagnetic configuration suggested in Ref. i45i, where all 
spins are oriented either parallel or antiparallel to the 
c direction. In both cases, the relative orientation of 
the four symmetry-related magnetic moments on the Fel 
sites or the Fe2 sites, respectively, is dictated by the mag- 
netic space group Pc'2[n. 

We note that in the magnetic configuration used by 
Arima et al. (Ref. lisl ) the net magnetization stems 
mainly from the intersite disorder; the system is a per- 
fectly compensated antiferromagnet if (i) all Fel and Fe2 
sites are occupied by Fe cations, (ii) there is no Fe occu- 
pation of the Gal and Ga2 sites, and (iii) the magnetic 
moments on the two Fe sites are the same. Here, we con- 
sider only this perfectly compensated configuration with 
no site disorder and composition a; = 1, i.e. all Fel and 
Fe2 sites are occupied by Fe^+ cations and all Gal and 
Ga2 sites are occupied by Ga'^"'' cations. Thus, the mag- 



TABLE III: Positions and magnetic moment directions rrii 
of the Fe sites in GaFeOa with magnetic space group Pc'2'in. 
Both Fe sites correspond to Wyckoff positions 4a. The mag- 
netic configuration is that discussed in Ref. [i^ . S is the mag- 
nitude of the Fe magnetic moment, which is assumed to be 
identical on both sites and a, b, c axe the usual orthorhombic 
lattice parameters. 
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netic configuration discussed in the following does not 
have a net magnetization. We point out that in gen- 
eral the magnetic and toroidal properties will depend on 
the exact occupation numbers of the various cation sites. 
The positions of the Fe sites as well as the corresponding 
moment directions are listed in Table [1111 

The magnetic space group Pc'2[n breaks both space 
and time reversal symmetries and thus allows the exis- 
tence of a macroscopic toroidal moment. Evaluation of 
Eqs. ([7|) and (|16p for the positions and moment directions 
listed in Table [1111 leads to the toroidization: 



2n 



^4(yi 
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-y2)b 

ma 




nb 



(24) 



It can be seen that there is a nontrivial toroidization 
along the a direction, as dictated by the magnetic space 
group symmetry, whereas the component along the b 
direction represents only the trivial increment resulting 
from the periodic boundary conditions, and the compo- 
nent along c is zero. The macroscopic toroidization along 
a depends on the difference of the coordinates yi and ?/2 
of the two different Fe sites along b. 

One can verify that for yi — y2 = jl (for any inte- 
ger I) the "magnetic lattice", i.e. the spatial arrange- 
ment of magnetic moments on the Fel and Fe2 sites, is 
centrosymmetric, and thus the system is non-toroidic in 
the localized moment limit (see also Ref. [ish . This is 
consistent with Eq. ([24]), which for yi — 2/2 = \i results 
only in a trivial non-toroidal component of T along the 
a direction. We point out that the non-toroidicity for 
yi — y2 = \l holds true only for the case of localized 
magnetic moments, where the presence of all nonmag- 
netic ions is neglected. The full crystallographic symme- 
try of this system (given by both magnetic and nonmag- 
netic ions) is non-centrosymmetric even for yi—y2 = \l- 
In fact, GaFeOs is an example of a pyroelectric crys- 
tal that is polar but not ferroelectric, i.e. the polariza- 
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tion cannot be switched, since it does not result from 
a small distortion of a centrosymmetric reference struc- 
ture. However, in the localized moment picture the sys- 
tem is non-toroidal for yi — y2 — jl and we can eval- 
uate the spontaneous toroidization with respect to this 
reference configuration. The relevant structural param- 
eters determined experimentally in Ref. at 4K are: 
a = 8.719 A, b = 9.368A, c = 5.067 A, yi = 0.5831, 
and j/2 = 0.7998. This gives a spontaneous toroidization 
of Ts = ^(0.25 - 0.2167) = 7.5 • 10-VB/A^ corre- 
sponding to a spontaneous toroidal moment per unit cell 
of 3.1/xbA. Here, we used the formal magnetic moment 
S = 5^iB of Fe3+. 

It can be seen from Eq. ([M]) that the set of toroidiza- 
tion values is also centrosymmetric around the origin if 
2/1 — 2/2 is equal to any integer multiple of | , even though 
only for yi — 2/2 — jl the corresponding magnetic mo- 
ment configuration is non-toroidal as discussed in the 
previous paragraph. This shows that even though the 
set of toroidization values of a non-toroidal structure is 
always centrosymmetric, the converse is not necessarily 
true. A centrosymmetric set of toroidization values does 
not necessarily correspond to a non-toroidal state. Thus, 
for 2/1—2/2 = (2fc -I- 1) • I the higher symmetry of the 
toroidization values is accidental and does not correspond 
to a vanishing macroscopic toroidization. 

The toroidal moment of a single unit cell of GaFeOa 
was also evaluated in Ref. [23l. without taking into account 
the multivaluedness of the toroidization due to the peri- 
odic boundary conditions. A value of to = 24.155 ^^A 
along the a direction was reported for the centrosym- 
metric reference structure with j/i — 2/2 = —0.25, and the 
spontaneous toroidal moment was specified as 0.03 io. It 
is unclear from Ref. [2^ whether mixing of Fe ions onto 
the Ga sites was included in the calculation and therefore 
a direct comparison with our calculation is not possible. 



D. BiFeOa 



bulk spiral component as well as the small weakly fer- 
romagnetic component from the discussion. Depending 
on the direction of the antiparallel Fe moments, the mag- 
netic point group is either 3to (if the moments are aligned 
along the polar z axis), m (if the moments are oriented 
perpendicular to the polar axis and parallel to the glide 
plane), or m' (with the moments perpendicular to the 
glide plane); all of these symmetries allow for a macro- 
scopic toroidization. First-principles calculations showed 
that 3m symmetry, which would not allow for a macro- 
scopic magnetization, is energetically unfavorablei^ 

BiFeOs provides an instructive example to illustrate 
the limitations of the localized moment approach to 
toroidal moments based on Eq. ([7]). Both space and 
time- inversion symmetries are broken in BiFeOs, so that 
in principle a toroidal moment is symmetry-allowed for 
this system. Indeed, an antisymmetric component of 
the magneto-electric tensor, indicating a non-vanishing 
toroidal moment, has been measured at high magnetic 
fields where the bulk spiral spin structure is destroyedi^ 
However, if we consider only the localized spins on the Fe 
sites, the magnetic lattice is centrosymmetric, i.e. it has 
a higher symmetry than the full magnetization density 
/x(r) (see the discussion towards the end of Sec. |TT|, and 
thus the corresponding toroidal moment vanishes. In the 
i?3c crystal structure of BiFeOs the inversion symme- 
try is lifted by displacements of the different ionic sub- 
lattices relative to each other. Since the local moment 
picture neglects the presence of all nonmagnetic ions, this 
inversion symmetry-breaking is not present in the purely 
magnetic lattice. 

If we take the rhombohedral axis of BiFeOa as the z- 
direction (hexagonal setup of the rhombohedral unit cell) 
and consider a perfect G-type ordering with the magnetic 
moments of the Fe cations along ±x, then the calculated 
toroidization is: 



BiFeOs is a multiferroic material of high practical in- 
terest since it combines both magnetic and ferroelectric 
order above room temperature (see Ref. l46h . It exhibits 
a rhombohedrally distorted perovskite structure (space 
group i?3c) involving both polar displacements of the ions 
along the [111] direction and counter- rotations of oxygen 
octahedra around this direction41i^i^ The spin struc- 
ture of BiFeOs is a superposition of various components. 
In a first approximation, the spins order in a G-type anti- 
ferromagnetic structure, where all neighboring magnetic 
moments are oriented antiparallel to each other In ad- 
dition, in bulk BiFeOs the axis along which the spins 
arc aligned rotates throughout the crystal, leading to an 
additional spiral spin structure with a large period of 
^620 A. However, this spiral component is absent in 
thin film sampleSf^ where instead a weak ferromagnetic 
moment, resulting from a small canting of the magnetic 
moments, has been report ed j'^^'"^^ Here, we exclude the 
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(25) 



Here, c is the lattice constant of BiFeOs along z within 
the hexagonal setup. The first term along the y direction 
in Eq. (j25|) is the trivial half-toroidization increment, and 
the whole set of toroidization values is centrosymmetric, 
i.e. non-toroidal. This shows that it is not always suf- 
ficient to consider magnetic dipole moments localized at 
the sites of the magnetic cations, but that the spatial 
moment distribution can be very important. The exact 
magnetic moment density /i.(r) always reflects the full 
magnetic space group symmetry of the system, whereas 
the reduction to localized magnetic moments can result 
in a higher symmetry than that of the full system. 
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VII. SUMMARY, CONCLUSIONS, AND 
OUTLOOK 

In summary we have presented a detailed study of 
magnetic toroidal moments in bulk periodic solids in the 
limit where the toroidal moment is caused by a time and 
space reversal symmetry breaking arrangement of local- 
ized magnetic moments. We have reviewed the basic mi- 
croscopic definitions and showed that the periodic bound- 
ary conditions lead to a multivaluedness of the toroidiza- 
tion, which suggests that only differences in toroidization 
are well-defined observable quantities. We suggest that 
the origin dependence of the toroidal moment should be 
treated by decomposing the magnetic moment arrange- 
ment into a fully compensated antiferromagnetic and an 
uncompensated ferromagnetic component, so that only 
the ferromagnetic component depends on the origin. Dif- 
ferences in toroidization resulting from the compensated 
part of the moment configuration can be evaluated rather 
straightforwardly, if one keeps in mind the macroscopic 
symmetry properties of the system. We have illustrated 
the main concepts and difhculties in evaluating magnetic 
toroidization in periodic systems by first discussing the 
simple example of a one-dimensional antiferromagnetic 
chain, and we have then analyzed the toroidization for 
four example materials. 

In addition to illustrating the general consequences of 
the origin dependence and the multivaluedness of the 
toroidal moment in periodic systems in the localized mo- 
ment limit, our main conclusion is that it is important to 
be aware of the macroscopic symmetry properties when 
evaluating toroidization changes. This is particularly 
striking in the example of the distorted one-dimensional 
antiferromagnetic chain discussed in Sec|Vl where the 
change in polarization due to a structural distortion can 
be calculated straightforwardly, whereas in the case of a 
magnetic moment reversal one has to subtract the im- 
proper toroidization change that is caused by the corre- 
sponding change in the toroidization increment. We have 
also shown the limitations of the local moment picture in 
evaluating toroidal moments in cases where the reduc- 
tion to localized moments changes the symmetry of the 
system. 

An open question, which we have not addressed in our 
theoretical analysis, is how the spontaneous toroidization 
can be measured experimentally. According to the funda- 
mental definition of the toroidal moment, this is in prin- 
ciple possible by measuring the torque on a sample that is 
placed in an inhomogeneous magnetic field. However, ei- 
ther a field with a constant curl over the whole dimension 
of the sample has to be generated, or the effects of other 
multipole moments that couple to other field components 
have to be subtracted appropriately. To our knowledge 
such a measurement has not been attempted as of yet. So 
far, experimental evidence for toroidal moments has been 
based mostly on the detection of an asymmetric compo- 
nent of the magneto-electric tensor, but since the corre- 
sponding prefactors are not known, an absolute quanti- 



tative determination of T is not possible. Similarly, the 
nonlinear optical techniques used in Ref. 7 are mostly 
sensitive to symmetry breaking, but are at best semi- 
quantitative. Overall it appears that quantitative mea- 
surements of toroidal moments are a challenging task, 
but in principle possible. 

Finally, we comment on some aspects of toroidal mo- 
ments that have been discussed in a sometimes confusing 
way in the literature. First, the relation between the 
toroidal moment and the asymmetric component of the 
linear magneto-electric effect, and second the outer prod- 
uct of the polarization with the magnetization, which has 
sometimes been interpreted as a toroidal moment. 

From a macroscopic symmetry point of view, the sym- 
metries which allow for a macroscopic toroidal moment 
are identical with that allowing for an antisymmetric 
component of the linear magneto-electric effect tensor. 
The relation between these two quantities can be seen by 
analyzing the following free energy expression (see also 
Ref.ini): 

U = -kP^ P E + -BM^ -M H 
2 2 

+ -AT^ + -CT^ + aT-(PxM) . (26) 
2 4 

This is the simplest possible free energy expression that 
can simultaneously describe (i) a phase transition from 
a para-toroidic (T = |T| = 0) into a ferrotoroidic phase 
(T 7^ 0), (ii) the coupHng of the electric polarization P 
and the magnetization M to the electric field E and the 
magnetic field H, respectively, and (iii) a coupling be- 
tween the electric polarization, the magnetization, and 
the toroidization. In Eq. ([26]) k and B are the inverse 
electric and magnetic susceptibilities, A and C are tem- 
perature dependent coefficients, and a determines the 
strength of the magneto-electric coupling. The trilinear 
form of the coupling term in Eq. ((26)) is the lowest pos- 
sible order that is compatible with the overall space and 
time reversal symmetries. The equilibrium values for P 
and M can be obtained by minimizing Eq. I|26p . This 
leads to: 

P = -{E - a{M xT)} (27) 

K 

and 

M = —{H -a{T X P)} . (28) 
B 

If one inserts Eq. ([^5]) into Eq. one obtains (to lead- 
ing order in T): 

P=-E--^{HxT) . (29) 

The last term in Eq. (j29[) represents an antisymmetric lin- 
ear magneto-electric effect proportional to the toroidiza- 
tion. Thus, the presence of the trilinear coupling term 
between toroidization, magnetization, and polarization 
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in Eq. (|26p gives rise to an antisymmetric magneto- 
electric effect P — aH in the ferrotoroidic phase, with 
'^v ~ ^ijkTk- Note that in general other terms 

in the free-energy expansion can give rise to additional 
symmetric contributions to the magneto-electric tensor, 
which are not proportional to the toroidal moment. 

In Eq. (I26p only the magnetization and the polariza- 
tion couple to H and E, the toroidization in general 
does not couple to any homogeneous external fields, in 
agreement with the fundamental definitions discussed in 
Sec. In] (in particular Eq. ([3])). The effective couphng 
represented by the invariant -Eme ^ T ■ [E x H) (dis- 
cussed in Refs. arises from the trilinear coupling 
term in Eq. (|26p if P and M are substituted by their 
corresponding fields, by inserting Eqs. ([27]) and (f28|) into 
Eq. pS)) . Of course it depends on the problem at hand 
whether it is more convenient to use the fields E and 
H or the magnetization M and polarization P as free 
variables. It is worth pointing out, though, that the two 
cases should be carefully distinguished. One can either 
use a description where P, T, and M are the free vari- 
ables, P and M couple linearly to their corresponding 
fields, and there is a trilinear coupling term of the form 
T ■ {P X M) , or one can alternatively use a picture where 
the variables P and M are eliminated altogether and are 
replaced by the field variables E and H. In this case the 
trilinear coupling term leads to an effective coupling of 
T to E X H. 

Another source of confusion is the relation of the 
toroidization to the cross product P x M, which has 
the same time and space reversal symmetry as T. This 
has led to a number of instances in the literature in 
which P X M itself has been described as a "toroidal 
moment".— We point out that this is generally not 
correct. From our discussion in Sec. Illll it becomes clear 
that P X M can never describe a toroidal moment re- 
sulting from the compensated part of a magnetic moment 
configuration (for which M — 0). We have also shown 
in Sec. IIIII that the change in toroidal moment resulting 
from the uncompensated part of the magnetic moment 
configuration can be expressed as At = 1/2 Afi x m, 
where AR is the average (non-centrosymmetric) dis- 
placement of the magnetic moments. We point out that, 
at least in the localized moment picture, AR is in general 
not proportional to P, and therefore the toroidal moment 
is in general not proportional to f x M . In most mag- 



netic ferroelectrics, the polarization is related to a rigid 
shift of the magnetic cations relative to the non-magnetic 
ions, which does not affect the average position R of the 
magnetic cations. It is therefore not clear whether the 
conditions for T ^ P x M are fulfilled in any currently 
known material. 

It has been argued that ferrotoroidicity is a key concept 
for fitting all forms of ferroic order in a simple fundamen- 
tal scheme based on the different transformation proper- 
ties of the corresponding order parameters with respect 
to time and space inversion (see Refs. 030, in particular 
Fig. 2 in Ref. 7). The four fundamental forms of ferroic 
order listed in these references are ferroelasticity, ferro- 
electricity, ferromagnctism, and ferrotoroidicity, with or- 
der parameters transforming according to the four differ- 
ent representations of the "parity group" , which is gen- 
erated by the two operations of time and space reversal]^ 
A similar scheme has also been proposed in Ref. [^, but 
with the electric toroidal moment g (see Sec. U and also 
Refs. [Hill) as the time and space symmetric order pa- 
rameter instead of the ferroelastic strain tensor. The lat- 
ter classification scheme seems more natural to us, since 
in this case all ferroic order parameters are vector quan- 
tities. It is thus very important to clearly distinguish be- 
tween magnetic ferrotoroidicity and electric ferrotoroidic- 
ity, which correspond to different representations of the 
parity group. On the other hand, the existence of fer- 
roelasticity, with the second-rank strain tensor as time 
and space invariant order parameter, raises the question, 
whether other ferroic second-rank tensor order param- 
eters, corresponding to different representations of the 
parity group, can be identified in the future. 
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